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Probing the existence of the ^ pea k-^i S0 correlation in long Gamma 
Ray Bursts 
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ABSTRACT 

We probe the existence of the -Epcak-^iso correlation in long GRBs using a sample of 442 
BATSE bursts with known E'pcak and with redshift estimated through the lag-luminosity cor- 
relation. This sample confirms that the rest frame peak energy is correlated with the isotropic 
equivalent energy. The distribution of the scatter of the points around the best fitting line is 
similar to that obtained with the 27 bursts with spectroscopic redshifts. We interpret the scatter 
in the -E pca k — £\so plane as due to the opening angle distribution of GRB jets. By assum- 
ing that the collimation corrected energy correlates with -E pca k we can derive the observed 
distribution of the jet opening angles, which turns out to be log-normal with a peak value of 
-6.5°. 
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1 INTRODUCTION 



Since the discovery of the cosmological distance scale to Gamma 
Ray Bursts in 1997 (Djorgovski et al. 1997) several GRB red- 
shifts have been spectroscopically measured 1 . The small sample 
of bursts, now comprising 42 events, allowed in the latest years the 
investigation of the GRB rest frame properties. Several correlations 
have been discovered between GRBs spectral properties: 

(i) the spectral lag - luminosity correlation r-Z/i so (Norris, 
Marani & Bonnel 2000; Norris 2002) 

(ii) the variability - luminosity correlation V-Li so (Fenimore & 
Ramirez-Ruiz 2000; Reichart et al. 2001). 

(iii) the peak energy - isotropic energy correlation i5 P eak-i?iso 
(Amati et al. 2002; Lloyd & Ramirez-Ruiz 2002). We will call it 
the Amati correlation. 

(iv) the peak energy - isotropic peak luminosity correlation 
-Epcak-iiso.poak (Yonetoku et al. 2004). 

(v) the peak energy - collimation corrected energy correlation 
-E P cak--B 7 (Ghirlanda, Ghisellini, Lazzati 2004, GGL04 hereafter). 
We will call it the Ghirlanda correlation. 

All these correlations can be represented by powerlaws: their 
slopes are summarized in Tab. 1. Note that, for the last three corre- 
lations, the peak energy -E poa k of the vF v spectrum is obtained by 
fitting the GRB spectrum time-integrated over the burst duration. 

Despite several attempts to explain these correlations (see e.g. 
Eichler D. & Lenvison A. 2004; Liang E., Dai Z. & Wu X.F. 2004), 
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1 http://www.mpe.mpg.de/~jcg/grbgen.html for a continuously updated 
table. 



there is no universally accepted interpretation yet. One of the ma- 
jor drawback of these correlations is that they were found for a 
small sample of bursts, mainly due to the limited number of reliable 
redshift measurements. However, one test which can be performed 
even without knowing the redshift is to check if GRBs with known 
fluence and -Bp Ca k are consistent with - say - the Amati correlation, 
by considering all possible redshifts. 

This kind of consistency check has been performed very re- 
cently by two groups which tested the Amati and the Ghirlanda 
correlation. Nakar & Piran (2004, hereafter NP04) tested the Amati 
correlation in the observed plane -Epcak,obs--F (where F is the burst 
observed fluence) using a sample of bright BATSE bursts (from 
Band et al. 1993; Jimenez, Band & Piran 2000). They concluded 
that ~40% of the BATSE long (bright) bursts are inconsistent with 
this relation. 

Band & Preece (2005, hereafter BP05) extended this work 
with a larger sample of BATSE GRBs (760 bursts from Mallozzi 
et al. 1998) and tested both the Amati and the Ghirlanda correla- 
tion. They claim that at least 88% of their bursts are inconsistent 
with the Amati correlation and 1.6% with the Ghirlanda correla- 
tion. They also claim that these conclusions are robust, since they 
are model-independent. 

An important point to notice is what is meant by Amati corre- 
lation. Originally, Amati et al. (2002) found -E^k oc with a 
small scatter, but using only 9 bursts. Later, GGL04 (see their Tab. 1 
and Tab. 2) collected 23 GRBs (including the above 9) with redshift 
measurements and published peak energy, finding a slightly differ- 
ent slope and, more importantly, a larger scatter around the cor- 
relation, whose significance was nevertheless improved (due to the 
increased number of points). In the same paper, it was demonstrated 
that for all bursts with a good measurement of the jet opening angle 
(through the detection of a break in the light curve of the afterglow), 
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^iso.peak 






Ref 


T 

V 

-^pcak 


o.87±g;ll 

3+0.11 
0.5 ±0.05 


0.52 ±0.06 


0.706 ±0.04 


(1), (2) 
(3) 

(4), (5), (6) 



Table 1. Slopes of the correlations between rest frame GRB quantities (for 
instance: r <x i? * 7 pcak ). Ref: (1) Norris et al. 2000; (2) This work (Sec.3) 
(3) Fenimore & Ramirez-Ruiz 2000, Reichart et al. 2001; (4) Yonetoku et 
al. 2004; (5) Amati et al. 2002; (6) Ghirlanda, Ghisellini & Lazzati 2004. 



the scatter of the Amati correlation was entirely due to the differ- 
ent opening angles of the jet: bursts with the same collimation- 
corrected energy E-, and the same intrinsic S pca k have different 
Eiso because they have different jet opening angles. The scatter 
of the Amati correlation, far from being due, i.e. to ill-measured 
quantities, has a physical origin and is revealing of the distribution 
of the opening angles of the jets. 

Bearing this in mind, in this letter we test if -E pca k correlates 
with i?i so even if we do expect a large scatter around the best fit- 
ting line (if any). To this aim we use a large sample of bursts with 
"pseudo redshifts" which have been estimated through the lag- 
luminosity relation (Sec. 2). Having found that indeed the S pca k- 
Ei so correlation exists (Sec. 3), we follow GGL04 and interpret the 
scatter as due to the distribution of jet opening angles, finding it 
(Sec. 4). Finally (Sec. 5), we discuss why NP04 and BP05, who 
also use a large sample of BATSE bursts, arrive to different conclu- 
sions. 

During the completion of this work, we received a manuscript 
from Bosnjak et al. By studying the same problem with an inde- 
pendent and complementary approach, they arrive at the same con- 
clusions presented in this letter. 

In the following we adopt a standard ACDM cosmology with 
n A = 0.3, Q, m = 0.7 and h = 0.7. 



2 THE GRB SAMPLE 

As described above the Amati correlation involves the peak energy 
of the time integrated spectrum of the prompt emission. This is ob- 
tained by fitting the GRB spectrum time integrated over the entire 
duration of the burst (see e.g. Band et al. 1993). 

Yonetoku et al. (2004; Y04 hereafter) have recently pub- 
lished a catalog of 689 long BATSE GRBs obtained through the 
spectral analysis of the high spectral resolution BATSE data. The 
main selection criterium adopted by Y04 is on the peak flux (i.e. 
> 2 x 10~ 7 erg cm -2 s" 1 ). Band, Norris & Bonnell (2004, BNB04 
hereafter) derived for the BATSE GRBs in the sample of Mallozzi 
et al. 1998 (the same used in BP05) their pseudo redshifts through 
the lag-luminosity correlation (Norris et al. 2000). This sample 
contains 1 165 bursts. We have selected 578 sources contained both 
in Y04 and in BNB04 to have a sample of bursts with a good mea- 
surement of i?p Ca k and an estimate of their redshift. For each se- 
lected burst, we took the fluence from the current BATSE GRB 
catalog. 2 We excluded from the final sample the sources with par- 
ticularly large uncertainties on the parameters which are used for 
the test that we present in the next sections: 

• 21 bursts with iJ poa k ^ 40 keV, since this value is very close 
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Figure 1. Redshift distribution of the sample of bursts with pseudo red- 
shift (filled histogram, 442 objects), compared to the redshift distribution 
of GRBs with measured spectroscopic redshift (red-hatched histogram, 27 
objects). Distributions are normalized. 



to the BATSE threshold, and this makes the fitting procedure highly 
uncertain (see e.g. Lloyd & Petrosian 2002); 

• 5 bursts with errors in £ pca k greater than _E pca k itself; 

• 6 bursts with errors on their fluence greater than their fluence. 
This sample contains 546 GRBs and it will be used for the test dis- 
cussed in Sec. 5. 

In addition, to the aim of testing the Amati correlation with a sam- 
ple of GRBs with reliable redshift estimates we excluded: 

• 94 bursts with errors in their estimated lag greater than the lag 
itself. 

• 10 bursts with redshift larger than 10. 

The final sample contains 442 GRBs and this is used in Sect.3. Al- 
though the estimated pseudo redshifts should be taken with care 
when considering specific objects, the distribution of redshifts for 
a large sample of bursts should be more reliable. In Fig. 1 we show 
that the distribution of the 442 pseudo redshifts of our sample is 
similar to the distribution of spectroscopically measured redshifts 
(27 objects). Even if not complete, this sample should be represen- 
tative of moderately bright BATSE bursts. 



3 TESTING THE EXISTENCE OF AN Speak - E ISO 
CORRELATION 

Amati et al. (2002) showed that the rest frame i5 pca k — E[ so are 
correlated in 9 long GRBs detected by BeppoSAX. They found 
a correlation coefficient of 0.949 with a chance probability P = 
5 x 10 -3 . This correlation was confirmed by GGL04 using a larger 
sample (objects listed in their Tab. 1 and 2) of 23 GRBs, with spec- 
troscopically measured redshifts and published peak energy. Here, 
we update this sample with 4 more GRBs whose redshifts have 
been recently measured (see Ghirlanda et al. 2005) 3 . With these 27 
bursts the Spearman rank correlation coefficient is r s — 0.87 and a 
chance probability P = 3.4 x 10 -9 . 

In Fig. 2 we show the rest frame -Ep ea k--Ej ao of these 27 GRBs 
(red symbols). The fit with a powerlaw model to these data points 
(weighting for the errors on both coordinates) is: 



2 http://cossc.gsfc.nasa.gov/batse/BATSE_Ctlg/index.html 



3 See also http://www.merate.mi.astro.it/~ghirla/deep/blink.htm 



3 



1000 



> 

0) 




100 



10 50 10 51 io 52 1 53 1 54 10 55 
E ls „ [erg] 

Figure 2. Rest frame peak energy -Ep Ca k versus isotropic energy E lso . Red 
filled circles are the 27 GRBs from Tab. 1 and 2 of GGL04 updated with 
the 4 new events with spectroscopic redshifts and the long-dashed red line 
is their best powerlaw fit. The black crosses are the 442 GRBs with redshift 
given by the lag-luminosity correlation (BNB04). The solid blue line rep- 
resents their best powerlaw fit weighted for the errors on both coordinates. 
The dot-dashed line is the fit to the complete sample (442 + 27 GRBs). 
The results of these fits are given in the text. Also shown (triple dot-dashed 
green line) is the Ghirlanda correlation. The shaded regions represent the la 
and 3cr width of the dispersion of the black crosses around their correlation, 
as represented by the blue solid line. 
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with a x 2 = 129.6/25 dof. The distribution of the scatter measured 
along the fitting line (i.e. the distances of the data points to the 
fitting line) of the 27 GRBs is shown in Fig. 3 by the red-hatched 
histogram. Fitting with a Gaussian we find a dispersion a — 0.22. 

For the 442 GRBs with pseudo redshifts, we have estimated 
the errors on the rest frame E[ so and i? pca k as follows. From 
Speak = -E pcak , obs (l + z) we derive 



peak 



E, 



pcak.obs 



+ 



1 + 2 



(2) 



Similarly from Ei BO — A-Kd 2 L F / {1 + z), where dz, is the luminosity 
distance and F is the GRB bolometric fiuence, we obtain: 



( S iso ) 



di 



1 + z 



+ 



(3) 



where d' L — ddL/dz. 

BNB04 do not give the error on the estimated redshifts. How- 
ever, it can be obtained from Li so = 47rd|_$, where <!> is the peak 
flux, as follows: 
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(4) 



The term ox iBo /Z/i so can be obtained from the lag-luminosity 
relation Li so = 10 s r n , using the error on the lag, i.e. a T given by 
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Figure 3. Solid histogram: Scatter of the 442 GRBs around their power- 
law fit (Eq. 6). The black solid line is the gaussian fit to this distribution. 
Red hatched histogram: scatter of the 27 GRBs with spectroscopic redshift 
around their powerlaw fit (Eq. 1). The red solid line is the gaussian fit. 



BNB04 in their table. However, also the error on the slope n and on 
the normalization B contribute to the error on the luminosity. Nor- 
ris et al. (2000) and BNB04 do not give the uncertainty on these 
parameters. For this reason we recomputed the lag-luminosity re- 
lation (with the data of Table 1 of Norris 2000) in the "barycenter" 
of the data points, namely tj, = 0.036 s and Li SO ,b = 3.36 x 10 52 
erg s _1 , in order to treat B and n as uncorrelated. 

We found exactly the same slope and normalization but we 
could now also find the uncertainties on these quantities, i.e. n ± 
cr n — 1.149±0. 205 and (in the coordinate system of the barycenter 
the normalization is null but its error is nonzero) i.e. B ± ob — 
0.0 ± 0.116. We can finally compute the error: 



\ -Zviso / 



= 2 In 10 



(^lo g (^))V|] + (^) (5) 



Fig. 2 shows the 442 GRBs with pseudo redshifts (black sym- 
bols) 4 . Interestingly, we find that these 442 GRBs still show a 
highly significant correlation between the rest frame E pea ^ and 
Eiso with a Spearman correlation coefficient r s = 0.7 and a chance 
probability P = 2.1 x 10~ 65 . The partial correlation coefficient, 
i.e. removing the effect of z, is r p — 0.65 and, therefore, indicates 
that the correlation is not determined by the common dependence 
of Eq. 2 and Eq. 3 on the term (1 + z). The fit (weighting errors on 
both coordinates) with a powerlaw to these data gives: 



E, 



peak 



100 keV 



(4.34 ± 0.05) 



Ei so 



7.5 x 10 52 erg 



(6) 



with a x 2 = 1735/440 dof. In this case the scatter distribution 
(filled histogram in Fig. 3) has a dispersion a — 0.22, consistent 
with the dispersion of the 27 GRBs with spectroscopic redshifts. 
It is remarkable that both the slopes of the two correlations (Eq. 

4 Note that i?i so is computed accounting for the band k-correction as de- 
scribed in Eq.2 of GGL04. 
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1 and Eq. 2 above) and their scatter distributions are very similar. 
However, the two correlations have different normalizations: the 27 
GRB with spectroscopic redshifts seem to define an envelope to the 
distribution of all GRBs of this sample. A possible reason for this 
behavior is that GRBs with spectroscopically measured redshifts 
tend to have fluences larger than the average, at least when compar- 
ing GRBs detected by BATSE or BeppoSAX. This fact has been 
pointed out also by NP04. 

We finally consider the total sample of bursts, i.e. the 27 GRBs 
with spectroscopic redshifts plus the 442 with pseudo redshifts, and 
find a Spearman correlation coefficient r s — 0.7 and P = 1.5 x 
10 -71 , and the partial correlation coefficient is r p = 0.66. The fit 
with a powerlaw gives: 

/ „ \ 0.51±0.01 



100 keV 



(3.64 ± 0.04) 



-Biso 

7.9 x 10 52 erg 



(V) 



with ax 2 = 2185/463 dof. 



4 THE ORIGIN OF THE SCATTER IN THE Speak-Siso 
PLANE 

Both the 27 GRBs with spectroscopic redshifts and the 442 GRBs 
here considered have a relatively large scatter in the E pBa k-Ei S o 
plane. In GGL04 it has been shown that for the few GRBs with 
measured 0j e t this scatter reduces to only 0.1 dex on the Ghirlanda 
correlation. This suggests (see also Firmani et al. 2005) that the 
scatter might be due to the angle distribution of GRBs. 

If we make the hypothesis that the Ghirlanda correlation has a 
scatter much smaller than the Amati correlation, we can derive the 
jet opening angle distribution for the 442 GRBs with pseudo red- 
shifts. This is compared to the angle distribution of the few GRBs 
with measured jet break time and redshift (see Tab. 2 of GGL04) in 
Fig. 4. The angle distribution derived for the 442 GRBs and repre- 
sented in Fig. 4 is described by a log-normal function: 



27T 0jetai 



exp 



(lnfljct - a 2 ) 2 
2a? 



(8) 



The best fit parameters are (ai, 02) = (2.2, 0.57) which lead to a 
peak value of ~6.5°. Note that the derived distribution is the sum 
of the distribution of angles for bursts with any value of S pC ak, 
and does not pretend to be representative of the distribution at each 
Speak- Indeed, if the jet opening angle distribution remains the 
same for all values of S pe ak, then we should expect the same slope 
for the Ghirlanda and Amati correlation. Since this is not the case, 
we can conclude that the current data indicate a trend: for smaller 
£" P eak and then E~ t , the average jet opening angle is larger than for 
large values of S pC ak and B 7 . It is premature to conclude whether 
this trend is due to selection effects or if it has a physical origin (i.e. 
smaller jet opening angles for intrinsically more powerful bursts), 
and we plan to investigate in more detail this problem in a forth- 
coming paper. 



5 COMPARISON WITH PREVIOUS WORKS 

In the previous section we have shown that if the pseudo-redshifts 
obtained through the lag-luminosity relation are indicative of real 
redshifts, then the -E P eak-Si so correlation exists even considering 
a large sample of GRBs. We now compare this conclusion with the 
conclusions of NP04 and BN05. 

The method adopted by NP04 and BP05 starts from a given 



0.1 



0.01 



0.001 




100 



iel [degrees] 



Figure 4. Jet opening angle distribution for the sample of GRBs with 
pseudo redshift (blue hatched histogram). This has been obtained requiring 
they obey the Ghirlanda correlation. The solid black line is the log-normal 
fit to this distribution, which peaks at ~ 6.5°. Also shown is the distribution 
(red hatched histogram) of measured jet opening angles for the 17 bursts in 
the sample of GGL04 of known redshifts and jet break time (Ghirlanda et 
al. 2005). 



rest frame correlation -E pca k = KE?^ where the two quantities 
are related to the observables Speak, obs = S pGa k/(l + z) and the 
fluence F = E iso (l + z) /And 2 L (z). We can form the ratio, 



E 



1/') 

P eak,obs 



47rd|(z) 



(1 + 2) 



(9) 



where a corresponds to the scatter of data points around the corre- 
lation that is being tested. The RHS of the above relation is a func- 
tion of z and 77 only, and is monotonically increasing with redshift 
if 77 ^ 0.75. For lower values of rj it has a maximum. 

As in the case of NP04 and BP05, we take into account now 
the upper limit of the ratio E^ B 2 k ohs /F. This upper limit estab- 
lishes an allowance region boundary on the corresponding plane. 
Furthermore, we include the scatter around the assumed correla- 
tion increasing the previous upper limit by the factor 10 CT ' 77 (see 
Eq. 9). 

NP04 and BP05 are testing the S pca k — Ei so correlation as 
originally found with 9 bursts by Amati et al. 2001 , i.e. with slope 
rj = 0.5. They allow for the possible scatter around this correlation 
with a factor 2 in the normalization of Eq. 9. This corresponds, for 
their assumptions, to consider a ~ 0.15. As shown in this work 
(and see also GGL04) the 27 GRBs show a dispersion with a — 
0.22 in the Speak — Siso plane around a correlation with slope 
rj — 0.56 (see Eq.l). Moreover, we consider the 3a consistency of 
our sample of bursts with this correlation. 

For the above reasons we applied the same test adopted by 
NP04 and BP05, but using the correlation found here, with its 
proper scatter, i.e a — 0.22. Note also that we use the peak en- 
ergy of the time integrated spectrum S pea k. 

In Fig. 5 we report the 27 GRBs with spectroscopic redshifts 
(open red circles, see. Sec. 3, plus 3 upper limits), as well as the 442 
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Figure 5. Observed plane. Fluence versus peak energy of the average spec- 
tram for the sample of 546 GRBs selected in Sec. 2. Black dots are the 442 
GRBs which has been used in Sec. 3. Grey stars are 104 GRBs excluded in 
the analysis of Sec. 3 either for their large redshifts or uncertain lags (see 
Sect.2). Also shown are the 30 GRBs (open red circles) with known red- 
shift (from GGL04, including/upper/lower limits on -E poa k). The dashed 
red line (black solid line) represents the boundary defined by the correlation 
of Eq. 1 (Eq. 6). The hatched red region (yellow filled region) represents 
the 3a shift of the corresponding boundary due to the scatter. Points above 
these boundaries are fully consistent with the corresponding correlations. 
Upper/lower limits for data points are indicated by arrows. The dot-dashed 
black line represents the boundary of the Amati correlation adopted in BP05 
and NP04. 



GRBs (black dots) of the sample described in Sec. 2. The dashed 
red line (black solid line) represents the boundary defined by Eq. 1 
(Eq. 6), while the red (yellow) hatched area shows the 3a shift of 
this boundary related to the scatter. For a gaussian distribution the 
area below the 3a limit is highly improbable for the elements of 
the corresponding sample. The point here is that 1.4% of the black 
dot distribution extends below the 3a limit of Eq. 1 correlation (see 
Tab. 2 for the la and 2a fractions). This result makes Eq. 1 in- 
compatible with the 442 BATSE GRB sample, in agreement with 
Nakar & Piran (2004) and Band & Preece (2005). This conclusion 
is strengthened by the fact that our allowance region is obtained 
using the upper limit of E^ k>oh jF. 

We conclude that Eq. 1 has to be considered as a temporary es- 
timate of the £'peak--Biso (Amati) correlation with an associated un- 
certainty likely due to the still small number of GRBs with known 
2. However, our main point is that the BATSE GRB sample con- 
firms the existence of an E pea ^-Ei so correlation and the scatter of 
this correlation (<j=0.22) agrees with the scatter of the sample of 
the 27 GRBs with spectroscopic redshifts (Eq. 1). 



6 CONCLUSIONS 

Through a large sample of 442 BATSE GRBs, with known spectral 
peak energy (Y04) and pseudo redshift (BNB04), we have com- 
puted the rest frame peak energy -E pca k and the k-corrected equiva- 
lent isotropic energy Ei so . We have found that these two quantities 
are correlated: the power-law fit to this correlation gives a slope 
slightly flatter and an i5i so a little smaller than the slope and the 
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2a 


3a 


Eq. 1 (red dashed) 


32% 


68% 


91.8% 


98.6% 


Eq. 6 (black solid) 


59% 


91.1% 


98.7% 


100% 



Table 2. Consistency of the sample of 546 GRBs described in Sec. 2 with 
the Amati correlations as given by Eq. 1 and Eq.6. The percentage represents 
the fraction of points above the lines in Fig. 5 (no scatter) or above the lower 
boundary of the corresponding a scatter. 

i?iao found with the 27 GRBs with spectroscopic redshifts, respec- 
tively. The scatter of the 442 data points around the fitting line is 
equal to the scatter of the 27 GRBs with known redshift, the exis- 
tence of outliers being marginal. 

We have further interpreted the scatter of these data points as 
due to the distribution of the jet opening angle of GRBs. By requir- 
ing that they obey the Ghirlanda correlation, we could estimate this 
distribution, which is a log-normal with an average angle of about 
6.5°. 

Our main conclusion is that the existence of an E pea ^-Ei so 
(Amati) correlation is confirmed, its scatter is probably rather well 
known, but its definitive estimate needs a wider sample of GRBs 
with spectroscopic measured redshifts. 
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